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Abstract 

These lectures are devoted to the low energy limit oi J\f = 2 SUSY gauge theories, which is 
described in terms of integrable systems. A special emphasis is on a duality that naturally acts on 
these integrable systems. The duality turns out to be an effective tool in constructing the double 
elliptic integrable system which describes the six-dimensional Seiberg-Witten theory. At the same 
time, it implies a series of relations between other Seiberg-Witten systems. 
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Part I 

Seiberg-Witten theory: gauge versus 
integrable theories 

The main target of these lectures is to present an integrable framework for one of the most intriguing 
discoveries of recent years in quantum field theory - a possibility of exact, non-perturbative treating 
the low energy behaviour of supersymmetric gauge theories proposed by N.Seiberg and E.Witten 
||. The results obtained turned out to be most effectively described by integrable systems Q. Since 
then, many new features of Seiberg-Witten theories have been discovered, all of them this or that way 
related to integrable systems (an intensive discussion of integrability in Seiberg-Witten theory and related 
structures (Whitham hierarchies and WDVV equations) can be found in the book Q that contains recent 
comprehensive reviews). Here we are going to concentrate on the most recent and puzzling structure, a 
duality present in integrable systems of the Toda/Calogero/Ruijsenaars family [||. 

Therefore, technically the lectures are divided into two parts. The first part is devoted to the Seiberg- 
Witten theories and their treatment within the integrable approach. The second part deals with the 
duality and its role in the Seiberg-Witten theories. 



1 Seiberg-Witten solution 

Being one of the most fascinating breakthroughs in field theory, the Seiberg-Witten approach opened 
new possibilities to deal with nonperturbative physics at strong coupling that amounted from the exact 
derivation of the low energy effective action for the TV = 2 SYM theory 0, In principle, to get the 
exact answer one needs to perform the summation of infinite instanton series. Therefore, the result by 
Seiberg and Witten was the first example of the exact derivation of the total instanton sum. However, 
the authors did not calculate instantonic series, using instead indirect and quite sophisticated methods. 
The calculation of the exact mass spectrum of stable particles was their second essentially new result. 
What is important, this theory is non-trivial in one loop and belongs to the class of asymptotically free 
theories. In fact, the situation is a little bit unusual, since, despite the fact that the direct instanton 
summation has not been performed yet, there are no doubts in the results derived which survives under 
different tests. 

In order to get the effective action unambiguously, the authors of [0, |^ used the three ideas: holomor- 
phy duality and their compatibility with the renormalization group flows. Holomorphy implies 
that the low energy effective action can be only of the form J Im.F(5') depending on a single holomorphic 
function J- called prepotential. The duality principle which provides the relation between the strong 
and week coupling regimes is natural in the UV finiteQ say JV = 4 SUSY |^ , theory, where the evident 
modular parameter built from the coupling constant and the O-teim 
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is not renormalized and coincides with its asymptotic value. In the theory with non-trivial renormaliza- 
tions, the duality transformation does no longer interchange the weak and strong coupling regimes of the 
same theory. Instead, it is much similar to the usual electric-magnetic duality and interchanges the low 
energy regimes of different theories. To construct this transformation explicitly, one needs an artificial 
trick, to introduce into the game an additional object, a higher genus Riemann surface whose period 
matrix coincides with the matrix of couplings depending on the values of condensate in a given vacuum. 

In the JV — 2 SYM theory this scheme is roughly realized as follows ||]. First of all, in the 
Lagrangian there is a potential term for the scalar fields of the form 

n0)=Tr[0,0+]2, (1.2) 

where the trace is taken over adjoint representation of the gauge group (say, SU{Nc)). This potential 
gives rise to the valleys in the theory, when [(/>, 0+] — 0. The vacuum energy vanishes along the valleys, 

^In fact, one should require more: not renormalized and, therefore, non-running coupling constants. 
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hence, the supersymmetry remains unbroken. One may always choose the v.e.v.'s of the scalar field to 
lie in the Cartan subalgebra (f) = diag(ai, an). These parameters can not serve, however, good 
order parameters, since there is still a residual Weyl symmetry which changes ai but leave the same 
vacuum state. Hence, one should consider the set of the gauge invariant order parameters Uk =< Tr0'^ > 
that fix the vacuum state unambiguously. Thus, we obtain a moduli space parametrized by the vacuum 
expectation values of the scalars, which is known as the Coulomb branch of the whole moduli space of 
the theory. The choice of the point on the Coulomb branch is equivalent to the choice of the vacuum 
state and simultaneously yields the scale which the coupling is frozen on. At the generic point of the 
moduli space, the theory becomes effectively abelian after the condensation of the scalar. 

As soon as the scalar field acquires the vacuum expectation value, the standard Higgs mechanism works 
and there emerge heavy gauge bosons at large values of the vacuum condensate. To derive the effective 
low energy action one has to sum up the loop corrections as well as the multiinstanton contributions. 
The loop corrections are trivial, since they are vanishing beyond the one loop due to the supersymmetry. 
The explicit summation over instantons has not been done yet, however, additional arguments allow one 
to define the effective action indirectly. 

The initial action of TV = 2 theory written in TV = 2 superfields has the simple structure S{'i') — 
Imr J Tr^*^ with t = ^ + ^ . Now the Af — 2 supersymmetry implies that the low energy effective action 
gets renormalized only by holomorphic contributions so that it is ultimately given by a single function 
known as prepotential 5*6// (vP) = Im J ^(^'). The prepotential is a holomorphic function of moduli Uk 
(or v.e.v.'s) except for possible singular points at the values of moduli where additional massless states 
can appear disturbing the low energy behaviour. 

Thus, the problem effectively reduces to the determination of one holomorphic function. If one 
manage to fix its behaviour nearby singularities, the function can be unambiguously restored. One of 
the singularities, corresponding to large values of v.e.v.'s, i.e. to the perturbative limit is under control 
(since the theory is asymptotically free). All other singularities are treated with the use of duality and of 
the non-renormalization theorems for the central charges of the SUSY algebra. A combination of these 
two ideas allows one to predict the spectrum of the stable BPS states which become massless in the deep 
non-perturbative region and are in charge of all other singularities. 

The duality transformation can be easily defined in the finite TV = 4 SUSY theory just as the modular 
transformations generated by t — > and r — > r + 1. This makes a strong hint that the duality can 
be related with a modular space of some Riemann surfaces, where the modular group acts. In fact, the 
naive application of duality meets some serious difficulties. The reason is that, in the asymptotically 
free theory, one has to match the duality with the renormalization group^. This is non-trivial, since now 
T depends on the scale which is supposedly involved into the duality transformation. The solution of 
this problem is that one is still able to connect the duality and modular transformations if considering 
different low energy theories connected by the duality. Then, the duality acts on the moduli space of 
vacua and this moduli space is associated with the moduli space of an auxiliary Riemann surface, where 
the modular transformations act. 

At the next step, one has to find out proper variables whose modular properties fit the field theory 
interpretation. These variables are the integrals^ of a meromorphic 1-form dS over the cycles on the 
Riemann surface, ai and af 



(where i, j = 1, TVc — 1 for the gauge group SU{Nc)). 

These integrals play the two-fold role in the Seiberg-Witten approach. First of all, one may calculate 
the prepotential and, therefore, the low energy effective action through the identification of and 

^In fact, there are more problems with duality even in Af = 2 SUSY theories. In particular, it unifies into one duality 
multiplet the monopole and gauge boson supermultiplets, while they have different spins, see, e.g., flll, b|. 
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dJ- /da with a defined as a function of moduli (values of condensate) by formula (1.3). Then, using the 
property of the differential dS that its variations w.r.t. moduli are holomorphic one may also calculate 
the matrix of coupling constants 



The second role of formula (^) is that, as was shown these integrals define the spectrum of the stable 
states in the theory which saturate the Bogomolny-Prasad-Sommerfeld (BPS) limit. For instance, the 
formula for the BPS spectrum in the SU (2) theory reads as 



Mn,m = \na{u) + maD{u)\^ (1.5) 

where the quantum numbers n, m correspond to the "electric" and "magnetic" states. The reader should 
not be confused with the spectrum derived from the low energy behaviour which fixes arbitrarily heavy 
BPS states. The point is that the BPS spectrum is related to the central charge of the extended SUSY 
algebra ^{m},{n} = ain^ + afm^ pO| and, therefore, has an anomaly origin. On the other hand, anomalies 
are not renormalized by the quantum corrections and can be evaluated in either of the UV and IR regions. 

Note that the column (a^, af) transforms under the action of the modular group SL{2, 2) as a section 
of the linear bundle. Its global behaviour, in particular, the structure of the singularities is uniquely 
determined by the monodromy data. As we discussed earlier, the duality transformation connects different 
singular points. In particular, it interchanges "electric" , ai and "magnetic" , af variables which describe 
the perturbative degrees of freedom at the strong and weak coupling regimes respectively. Manifest 
calculations with the Riemann surface allow one to analyze the monodromy properties of dual variables 
when moving in the space of the order parameters. For instance, in the simplest SU{2) case, on the 
u-plane of the single order parameter there are three singular points, and the magnetic and electric 
variables mix when encircling these points. Physically, in the theory with non-vanishing O-term (not pure 
imaginary r), the monopole acquires the electric charge, while the polarization of the instanton medium 
yields the induced dyons. 



2 Seiberg-Witten theory and theory of prepotential 

Although the both auxiliary objects, the Riemann surface and the differential dS have come artificially, 
an attempt to recognize them in SUSY gauge theories results into discovery of the integrable structures 
responsible for the Seiberg-Witten solutions. 

Now let us briefly formulate the structures underlying Seiberg-Witten theory. Later on, we often refer 
to Seiberg-Witten theory as to the following set of data: 

• Riemann surface C 

• moduli space Ai (of the curves C) , the moduli space of vacua of the gauge theory 

• meromorphic 1-form dS on C 

How it was pointed out in Q , this input can be naturally described in the framework of some underlying 
integrable system. 

To this end, first, we introduce bare spectral curve E that is torus y'^ = x'^ + g2X^ + gs for the UV-finite 
gauge theories with the associated holomorphic 1-form duj = dx/y. This bare spectral curve degener- 
ates into the double-punctured sphere (annulus) for the asymptotically free theories (where dimensional 
transmutation occurs): x ^ w + l/w, y ^ w — l/w, duj = dw/w. On this bare curve, there are given 
a matrix- valued Lax operator L{x,y). The corresponding dressed spectral curve C is defined from the 
formula det(L — A) = 0. 

This spectral curve is a ramified covering of E given by the equation 

P(A;a;,y)=0 (2.6) 
In the case of the gauge group G = SU{Nc), the function T' is a polynomial of degree Nc in A. 
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Thus, we have the spectral curve C, the moduh space A4 of the spectral curve being given just by 
coefficients of V. The third important ingredient of the construction is the generating l-form dS = Xdoj 
meromorphic on C ("^" denotes the equality modulo total derivatives). From the point of view of 
the integrable system, it is just the shortened action "pdq" along the non-contractible contours on the 



Hamiltonian tori. This means that the variables Ui in (1.3) are nothing but the action variables in 
the integrable system. The defining property of dS is that its derivatives with respect to the moduli 
(ramification points) are holomorphic differentials on the spectral curve. This, in particular, means that 

= du^ (2.7) 

where duji are the canonical holomorphic differentials^ Integrating this formula over _B-cycles and using 
that qd = dJ-/da, one immediately obtains ( |l.4[ ). 

So far we reckoned without matter hypermultiplets. In order to include them, one just needs to 
consider the surface C with punctures. Then, the masses are proportional to the residues of dS at the 
punctures, and the moduli space has to be extended to include these mass moduli. All other formulas 
remain in essence the same (see |0 for more details). 

The prepotential J- and other "physical" quantities are defined in terms of the cohomology class of 



dS, formula (|f .3|). Note that formula (1.4) allows one to identify the prepotential with logarithm of the 



r-function of the Whitham hierarchy ||12|: T = logr. 



3 Seiberg-Witten theory and integrable systems 

Thus, we have fixed the set of data, that is, a Riemann surface, the corresponding moduli space and 
the differential dS that gives the physical quantities, in particular, prepotential, and explained how it 
can be associated with an integrable system. In the case of the SU{2) pure gauge theory, which has 
been worked out in detail in [|l| this correspondence can be explicitly checked |^. In particular, one may 
manifestly check formulas (1.3), since all necessary ingredients are obtained in Q|. In this way, one can 
actually prove that this Seiberg-Witten theory corresponds to the integrable system (periodic Toda chain 
with two particles, see below). However, to repeat the whole Seiberg-Witten procedure for theories with 
more vacuum moduli is technically quite tedious if possible at all. Meanwhile, only in this way one could 
unquestionably prove the correspondence in other cases. 

In real situation, the identification between the Seiberg-Witten theory and the integrable system 
comes via comparing the three characteristics 

• number of the vacuum moduli and external parameters 

• perturbative prepotentials 

• deformations of the two theories 

The number of vacuum moduli (i.e. the number of scalar fields that may have non-zero v.e.v.'s) on the 
physical side should be compared with dimension of the moduli space of the spectral curves in integrable 
systems, while the external parameters in gauge theories (bare coupling constant, hypermultiplct masses) 
should be also some external parameters in integrable models (coupling constants, values of Casimir 
functions for spin chains etc). 

As for second item, we already pointed out the distinguished role of the prepotential in Seiberg- 
Witten theory (which celebrates a lot of essential properties, see e.g. jl^, ^). In the prepotential, the 
contributions of particles and solitons/monopoles (dyons) sharing the same mass scale, are still distin- 
guishable, because of different dependencies on the bare coupHng constant, i.e. on the modulus r of the 
bare coordinate elliptic curve (in the UV-finite case) or on the Aqcd parameter (emerging after dimen- 
sional transmutation in UV-infinite cases). In the limit r ioo {Aqcd — > 0), the solitons/monopoles 
do not contribute and the prepotential reduces to the "perturbative" one, describing the spectrum of 
non- interacting particles. It is immediately given by the SUSY Coleman- Weinberg formula pl| : 

*I.e. satisfying the conditions 
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reps RA 

Seiberg-Witten theory (actually, the identification of appropriate integrable system) can be used to con- 
struct the non-perturbative prepotential, describing the mass spectrum of all the "light" (non-stringy) 
excitations (including solitons/monopoles). Switching on Whitham times presumably allows one to 
extract some correlation functions in the "light" sector. 



The perturbative prepotential (3.8) has to be, certainly, added with the classical part. 



^class = -Tfl^ (3.9) 

since the bare action is S{^) = Imr J Tr^^. Here r is the value of the coupling constant in the deeply UV 
region. It plays the role of modular parameter of the bare torus in the UV finite theories. On the other 
hand, in asymptotically free theories, where the perturbative correction behaves like a?\og{a / Kqc d) , t 
is just contained in Aqcd, since it can be invariantly defined neither in integrable, nor in gauge theories. 

In fact, the problem of calculation of the prepotential in physical theory is simple only at the pertur- 
bative level, where it is given just by the leading contribution, since the /3-function in A/" = 2 theories 
is non-trivial only in one loop. However, the calculation of all higher (instantonic) corrections in the 
gauge theory can be hardly done at the moment]^ Therefore, the standard way of doing is to make 
an identification of the Seiberg-Witten theory and an integrable system and then to rely on integrable 
calculations. This is why establishing the "gauge theories ^ integrable theories" correspondence is of 
clear practical (apart from theoretical) importance. 

Now let us come to the third item in the above correspondence, namely, how one can extend the 
original Seiberg-Witten theory. There are basically three different ways to extend the original Seiberg- 
Witten theory . 

First of all, one may consider other gauge groups, from other simple classical groups to those being a 
product of several simple factors. The other possibility is to add some matter hypermultiplets in different 
representations. The two main cases here are the matter in fundamental or adjoint representations. 
At last, the third possible direction to deform Seiberg-Witten theory is to consider 5- or 6-dimensional 
theories, compactified respectively onto the circle of radius R5 or torus with modulus R^/Rq (if the 
number of dimensions exceeds 6, the gravity becomes obligatory coupled to the gauge theory). 

Now we present a table of known relations between gauge theories and integrable systems, see Fig.l. 



Let us briefly describe different cells of the table. 

The original Seiberg-Witten model, which is the Ad pure gauge SU{N) theory (in fact, in their papers 
[0, ^, the authors considered the SU{2) case only, but the generalization made in is quite immediate), 
is the upper left square of the table. The remaining part of the table contains possible deformations. Here 
only two of the three possible ways to deform the original Seiberg-Witten model are shown. Otherwise, 
the table would be three-dimensional. In fact, the third direction related to changing the gauge group, 
although being of an interest is slightly out of the main line. 

One direction in the table corresponds to matter hypermultiplets added. The most interesting is 
to add matter in adjoint or fundamental representations, although attempts to add antisymmetric and 
symmetric matter hypermultiplets were also done (see 16 for the construction of the curves and for 



the corresponding integrable systems). Adding matter in other representations in the basic SU{N) case 
leads to non-asymptotically free theories. 



Columns: Matter in adjoint vs. fundamental representations of the gauge group. Matter in 
adjoint representation can be described in terms of a larger pure SYM model, either with higher SUSY or 
in higher dimensional space-time. Thus models with adjoint matter form a hierarchy, naturally associated 
with the hierarchy of integrable models Toda chain ^ Calogero ^ Ruijsenaars ^ Dell ^ |l^, ^ 
^ 22 , I23I] . Similarly, the models with fundamental matter are related to the hierarchy of spin chains 



originated from the Toda chain: Toda chain ^ XXX ^ XXZ ^ XYZ [gj. 

Note that, while coordinates in integrable systems describing pure gauge theories and those with 
fundamental matter, live on the cylinder (i.e. the dependence on coordinates is trigonometric), the 

^In order to check the very ideology that the integrable systems lead to the correct answers, there were calculated first 
several corrections | ^^ . The results proved to exactly coincide with the predictions obtained within the integrable approach. 
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Figure 1: SUSY gauge theories <J=^ integrable systems correspondence. The perturbative limit is marked 
by the italic font (in parenthesis). 

coordinates in the Calogero system (adjoint matter added) live on a torus^. However, when one takes 
the perturbative limit, the coordinate dependence becomes trigonometric. 

Lines: Gauge theories in different dimensions. Integrable systems relevant for the description 
of vacua of d = 4 and d — ^ models are respectively the Calogero and Ruijsenaars ones (which possess 
the ordinary Toda chain and "relativistic Toda chain" as Inozemtsev's limits [p5|), while d = 6 theories 
are described by the double elliptic (Dell) systems. When we go from Ad (Toda, XXX , Calogero) 
theories to (compactified onto circle) 5d (relativistic Toda, XXZ , Ruijsenaars) theories the momentum- 
dependence of the Hamiltonians becomes trigonometric (the momenta live on the compactification circle) 
instead of rational. Similarly, the (compactified onto torus) Qd theories give rise to an elliptic momentum- 
dependence of the Hamiltonians, with momenta living on the compactification torus. Since adding the 
adjoint hypermultiplet elliptizes the coordinate dependence, the integrable system corresponding to 6d 
theory with adjoint^ matter celebrates both the coordinate- and momentum-dependencies elliptic. A 
candidate for "the elliptic Toda chain" was proposed in |p6| . 

Further in this part we show how the Riemann surfaces and the meromorphic differentials arise when 
considering a concrete integrable many-body system. We are going to explain in more details the simplest 
top left corner of the table and make some comments on the adjoint matter case in different dimensions. 
Then, after discussing duality, we shall come to the most general adjoint matter system, the right bottom 
cell in the table which is our ultimate target in these lectures. It can be, however, done only with the use 
of the coordinate-momentum dualityrl 



4 4(i pure gauge theory: Toda chain 

Thus, we start with the simplest case of the Ad pure gauge theory, which was studied in details in the 
original paper for the SU[2) gauge group and straightforwardly generalized to the SU[Nc) case in 

This system is described by the periodical Toda chain with period Nc H, ^ ^ whose equations of 
motion read 

^Since these theories are UV finite, they depend on an additional (UV-regularizing) parameter, which is exactly the 
modulus T of the torus. 

^Let us point out that by adding the adjoint matter ^a^always mean soft breaking of higher supersymmetries. 
*This duality has nothing to do with the duality of pi, H which we discussed earlier. 
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(4.10) 



% _ ^ 

dt at 

with periodic boundary conditions Pi+N^ — Pi^ Qi+N^ — Qi imposed. In physical system, there are Nc 
moduH and there are no external parameters {Aqcd can be easily removed by rescaling). In the Toda 
system, there are exactly Nc conservation laws. These conservation laws can be constructed from the 
Lax operator defined for any integrable system. The eigenvalues of the Lax operator do not evolve, thus, 
any function of the eigenvalues is an integral of motion. 

There are two different Lax representation describing the periodic Toda chain. In the first one, the 
Lax operator is represented by the Nc x Nc matrix depending on dynamical variables 



( Pi 

g5("32-gi) 





V 



w 



,i(«2-9i) 

P2 e 

,5(93-92) 



The characteristic equation for the Lax matrix 





5(93-92) 

-P3 






PN^ 



(4.11) 



r{X,w) = detAT^xA'. {C^^H - A) 
2;enerates the conservation laws and determines the spectral curve 



= 



(4.12) 



w + - = 2Pn^{X) 
w 

where Pn^ (A) is a polynomial of degree Nc whose coefficients are integrals of motion. 



(4.13) 



The spectral curve (4.13) is exactly the Riemann surface introduced in the context of SU{Nc) gauge 
theory. The integrals of motion parametrize the moduli space of the complex structures of the hyperelliptic 
surfaces of genus Nc — 1, which is the moduli space of vacua in physical theory. 

If one considers the case of two particles {SU{2) gauge theory), P{\) ~ — u, u ~ p^ — coshg, 
p = Pi = —p2, q = <li — 02- Thus we see that the order parameter of the SUSY theory plays the 
role of Hamiltonian in integrable system. In the perturbative regime of the gauge theory, one of the 
exponentials in coshg vanishes, and one obtains the non-periodic Toda chain. In the equation ( [1.13| ) the 
perturbative limit implies vanishing the second term in the l.h.s., i.e. the spectral curve becomes the 
sphere w = 2Pn^{X). 

After having constructed the Riemann surface and the moduli space describing the 4(i pure gauge 
theory, we turn to the third crucial ingredient of Seiberg-Witten theory that comes from integrable 
systems, the generating differential dS. The general construction was explained in the Introduction, this 
differential is in essence the "shorten" action pdq. Indeed, in order to construct action variables, one 
needs to integrate the differential dS — J^iPi'^Qi over Nc — 1 non-contractable cycles in the Liouville 
torus which is nothing but the level submanifold of the phase space, i.e. the submanifold defined by 
values of all Nc — 1 integrals of motion fixed. On this submanifold, the momenta pi are functions of the 
coordinates, qi. The Liouville torus in Seiberg-Witten theory is just the Jacobian corresponding to the 
spectral curve ( 4.13 ) (or its factor over a subgroup) . 

In the general case of a (^-parameter family of complex curves (Riemann surfaces) of genus g, the 
Seiberg-Witten differential dS is characterised by the property SdS = Y^^^-yduidvi, where dvi{z) are g 
holomorphic 1-differentials on the curves (on the fibers), while 5ui are variations of g moduli (along the 
base) . In the associated integrable system, Ui are integrals of motion and tt^ , some g points on the curve 
are momenta. The symplectic structure is 



dai A dp'l"''^ = dui A dvi{'Kk) 



(4.14) 



4=1 



The vector of the angle variables. 
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TTfe 



/ du>^ (4.15) 

k=l 

is a point of the Jacobian, and the Jacobi map identifies this with the g-th. power of the curve, Jac = C®^. 
Here duJi are canonical holomorphic difi^erentials, dvi — X]j=i '^^j §a ^^i' ^ome details on the symplectic 
form on the finite-gap solutions can be found in p9[ |. 

Technical calculation is, however, quite tedious. It is simple only in the 2-particle {SU{2)) case, when 
the Jacobian coincides with the curve itself. In this case, the spectral curve is 

iu + -^2(X^^u) (4.16) 
w 

while u ~ — coshg. Therefore, one can write for the action variable 

dw .„ , dw 



pdq ~ i> \/u — cosh qdq — (f) A , dS — A — (4-17) 

' ' w w 



where we made the change of variable w = e'^ and used equation ( 4.16| ) 



Now one can naturally assume that this expression for the differential dS is suitable for generic 
Nc- A long calculation (which can be borrowed, say, from the book by M.Toda [0) shows that this 
is really the case. One can easily check that the derivatives of dS w.r.t. to moduli are holomorphic, 
up to total derivatives. Say, if one parametrize^ P/v^(A) = — A^'= + SAr^-2A^'=~^ + ... and note that 
dS = Xdw/w — XdPN^[X)/Y — Pf^^{X)dX/Y + total derivatives , then 

9d5 ^ A^ (4.18) 

dsk Y 

and these differentials are holomorphic if fc < Nc — 2. Thus, Nc — l moduli gives rise to iVc — 1 holomorphic 
differentials which perfectly fits the genus of the curve (we use here that there is no modulus sat^-i). 



5 Perturbative prepotentials 

In the previous section, we established in the concrete example the connection between the Seiberg-Witten 
integrable theory and the gauge theory. However, so far we studied only a kind of general landscape in the 
theory, since the real, numerical results enter the game only with the prepotential calculated. Moreover, 
the correspondence between gauge and integrable theories is supported by comparing the perturbative 
prepotentials. We discuss it now. 

The technical tool that allows one to proceed with effective perturbative expansion of the prepotential 
is "the residue formula" , the variation of the period matrix, i.e. the third derivatives of J- (a) (see, e.g. 
[0): 

d^J- dTi4 dojidujjduJk ,^ 

-reSd5=o jj^ , (5.1) 



daidajdak dak 6dS 

where SdS = d (^^) or, exphcitly 6dS = dXd(,. Here the differentia^^ d^ given on the bare torus 

degenerates into dw/w when torus degenerates into the punctured sphere (Toda chain). We remark that 
although d^ does not have zeroes on the bare spectral curve when it is a torus or doubly punctured sphere, 
it does in general however possess them on the covering C. 

In order to construct the perturbative prepotentials, one merely can note that, at the leading order, 
the Riemann surface (spectral curve of the integrable system) becomes rational. Therefore, the residue 
formula allows one to obtain immediately the third derivatives of the prepotential as simple residues on 



sphere. In this way, one can check that the prepotential has actually the form ( p. 81 ) |1 1| . Later we discuss 



these calculations, but right now we shall see what the perturbative prepotentials are in in gauge theory. 

^The absence of the term A^<=~^ is due to the SU{Nc) group and corresponds to the total momentum equal to zero, i.e. 
to the centre mass frame. This is important for the counting of genus. 

We reserve notation g for the variable living on the bare spectral curve (torus), while shall denote later through f the 
variable living on the compactification torus in 6 dimensions. 
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As a concrete example, let us consider the SU (Nc) gauge group. Then, say, the perturbative prepo- 
tential for the pure gauge theory acquires the form 



pert 



87ri 



)) log (a^-aj) 



(5.2) 



This formula establishes that when v.e.v.'s of the scalar fields in the gauge supermultiplet arc non- 
vanishing (perturbatively are eigenvalues of the vacuum expectation matrix (^)), the fields in the 
gauge multiplet acquire masses m^r' = — a^' (the pair of indices (r, r') label a field in the adjoint 
representation of SU{Nc))- The eigenvalues are subject to the condition '^i — 0- Analogous formula 
for the adjoint matter contribution to the prepotential is 



(5.3) 



These formulas can be almost immediately extended to the "relativistic" hd N = 2 SUSY gauge 
models with one compactified dimension. One can understand the reason for this "relativization" in the 
following way. Considering four plus one compact dimensional theory one should take into account the 
contribution of all Kaluza-Klein modes to each 4-dimensional field. Roughly speaking it leads to the 



1-loop contributions to the effective charge of the form (a^ 



T- 



E 

masses 



log masses ' 



E 



log a. 



m 



log (i?5a.y + to) log sinh i?5ai 



(5.4) 



i.e. coming from 4d to hd one should make a substitution aij smhR^aij, at least, in the formulas for 
perturbative prepotentials. Adding the adjoint matter in this case would give the effective charge (or, 
equivalently, the period matrix of the spectral curve of the Ruijsenaars-Schneider system |2^, see table 
1): 



log 



/ \/sinh {R^aij + e) sinh {R^ai 



sinh R^aij 



(5.5) 



The same general argument can be equally applied to the Gd case, or to the theory with two extra 
compactified dimensions, of radii i?5 and Rq. Indeed, the account of the Kaluza-Klein modes allows one 
to predict the perturbative form of charges in the Qd case as well. Namely, one should expect them to 
have the forrrf^ 



- E 

masses 

log I i?5a, 



masses ' 



E l m n 



] - log6'* 



R^aij 



.R, 
R( 



(5.6) 



i.e. coming from 4d {5d) to 6d one should replace the rational (trigonometric) expressions by the elliptic 
ones, at least, in the formulas for perturbative prepotentials. In this case adding the adjoint matter would 
lead to the effective charge 



log 



0* I R^Qij + e 



:f ) 9^[R,a,,-e 



Re 



6, (R^a,j 



(5.7) 



6 4(i theory: a simple example of calculation 

Here we describe the method of calculating prepotentials in the simplest example of the Toda chain. In 
the perturbative limit (""^^^ ~^ 0) the second term in ( [4.13D vanishes, the curve acquires the form 

^-"^We denote through 6* the S-function with odd characteristics [pl| . 
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W^2PnS^). PNA^)^X{{^-^i) (6.8) 



a rational curve with punctures |^ and the generating differential ( 4.17 ) turns into 



d4tlt = AdlogPjv.(A) (6.9) 
Now the set of the TVc — 1 independent canonical holomorphic differentials is 

■^•-(A^-A^)''^- (A-A'j(A''-A„J - ^-'■■■^•'V.-l. ^«-^.-^^ 

and the ^-periods ai = Xi (the independent ones are, say, with i — I, . . . , Nc — 1) coincide with the roots 
of polynomial Pn^{X), Xi. The prepotential can be now computed via the residue formula 

duJidujjdhJk 

^pert,zjk = res — 7TT-JT (6-11) 

which can be explicitly done. The only technical trick is that it is easier to compute the residues at the 
poles of du's instead of the zeroes of dlogP. This can be done immediately since there are no contributions 
from the infinity A = oo. The final results have the following form 



_^J_ _6_ ^ 1 

ZTT'lJ' pert,iii — / T ^ "7 ^" / "7 ^ 

k^'i k=tNc 

r> -T- 3 2 X - 1 XjJ\J 

2TnJ^p^rt,tzj ^ T + T + 2^ T T , » / J, (6.12) 



k^'i,jM, 



pert,ij 



giving rise to the prepotential formula 

^P-* = ^E/^'^(^^^)' = ^'log^' (6.13) 



This formula perfectly matches the perturbative prepotential expected in the gauge theory (5.2) 



Part II 

Dualities in the many-body systems and 
gauge theories 

In this part, we are going to discuss a property of some integrable systems which is called duality. This 
property will allow us to fill in the right bottom corner of table 1, i.e. to construct a double elliptic 
system. 

In fact, as in string theory there are several dualities (S-, T- and U-dualities) , as in integrable systems 
there are several kinds of duality [Q. Here we concentrate on the most interesting, important and 
still remaining mysterious one, the coordinate-momentum duality. This is a duality between pairs of 
dynamical systems (which is self-duality for some systems) ||, |3^, 34, |3|, |2^. Note that this duality 
at the moment is attributed only with the Hitchin type systems, i.e. the Seiberg-Witten theories with 
adjoint matter, while spin chains are still not embedded into this treatment. 

Let us consider a system with the phase space being a hyperkahler four dimensional manifold. Let this 
manifold involve at most two tori, i.e. be the elliptically fibered K3 manifold. One torus is that where 

^^These punctures emerge as a degeneration of the handles of the hyperelhptic surface so that the A-cycles encircle the 
punctures. 
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the momenta live and the other one is that where the coordinates hve. The duahty at hands actuaUy 
interchanges momentum and coordinate tori and, therefore, can be called pg-duality. Note that typically 
one deals with degenerations of these tori. 

One can also consider a theory which is self-dual. This integrable system is called Dell (double elliptic) 
system^ and describes the Seiberg-Witten theory for the six-dimensional theory with adjoint matter. 

The other cases correspond to some degenerations of Dell system. Say, the degeneration of the 
momentum torus to C jZ^ leads to the five-dimcnsional theory, while the degeneration to leads to the 
four-dimensional theory. Since the modulus of the coordinate torus has meaning of the complexified bare 
coupling constant in field theory, interpretation of the degenerations of the coordinate torus is different. 
In particular, the degeneration to the cylinder corresponds to switching off the instanton effects, i.e. to 
the perturbative regime. 

A generic description of the pq-duality expresses a relationship between two completely integrable 
systems 5'i, 5*2 on a fixed symplectic manifold with given symplectic structure (M, a;) and goes back to 
. We say the Hamiltonian systems are dual when the conserved quantities of 5*1 and 5*2 together form 
a coordinate system for M. Consider for example free particles, i?^^' = '^iP\/k. For this system the 
free particles momenta are identical to the conserved quantities or action variables. Now consider the 
Hamiltonian H^^^ — '^^Qi/k with conserved quantities g^. Together {pi,qi} form a coordinate system 
for phase space, and so the two sets of Hamiltonians are dual. Duality then in this simplest example 
is a transformation which interchanges momenta and coordinates. For more complicated interacting 
integrable systems finding dual Hamiltonians is a non-trivial exercise. Note that this whole construction 
manifestly depends on the particular choice of conserved quantities. A clever choice may result in the 
dual system arising by simply interchanging the momentum and coordinate dependence, as in the free 
system. 

Some years ago Ruijsenaars 1^ observed such dualities between various members of the Calogero- 
Moser and Ruijsenaars families: the rational Calogero and trigonometric Ruijsenaars models were dual 
to themselves while trigonometric Calogero model was dual with the rational Ruijsenaars system (see 
^ for more examples). These dualities were shown by starting with a Lax pair L = L{p, q) and an 
auxiliary diagonal matrix A — A{q). When L was diagonalized the matrix A became the Lax matrix for 
the dual Hamiltonian, while L was a function of the coordinates of the dual system. Dual systems for a 
model possessing a Lax representation are then related to the eigenvalue motion of the Lax matrix. We 
discuss this approach in more details below. 

Another approach to finding a dual system is to make a canonical transformation which substitutes 
the original set of Poisson-commuting coordinates q^, {g^, Qj} = 0, by another obvious set of the Poisson- 
commuting variables: the Hamiltonians hi{p,q) or, better, the action variables ai{h) = ai{p,q). It will 
be clear below that in practice really interesting transformations are a little more sophisticated: hi 
are identified with certain functions of the new coordinates (these functions determine the Ruijsenaars 
matrix A{q)), which - in the most interesting cases - are just the same Hamiltonians with the interactions 
switched-off. Such free Hamiltonians are functions of momenta alone, and the dual coordinates substitute 
these momenta, just as one had for the system of free particles. 

The most interesting question for our purposes is: what are the duals of the elliptic Calogero and 
Ruijsenaars systems ? Since the elliptic Calogero (Ruijsenaars) is rational (trigonometric) in momenta 
and elliptic in the coordinates, the dual will be elliptic in momenta and rational (trigonometric) in coordi- 
nates. Having found such a model the final elliptization of the coordinate dependence is straightforward, 
providing us with the wanted double-elliptic systems. 



7 j9g-duality: two-body system 

The way of constructing dual systems via a canonical transformation is technically quite tedious for 
many-particle systems, but perfectly fits the case of SU{2) which, in the center-of-mass frame, has 
only one coordinate and one momentum. In this case the duality transformation can be described 
explicitly since the equations of motion can be integrated in a straightforward way. Technically, given 
two Hamiltonian systems, one with the momentum p, coordinate q and Hamiltonian h(p, q) and another 
with the momentum P, coordinate Q and Hamiltonian H{P, Q) we may describe duality by the relation 

^''In fact, this is an example of the Dell system with 2 degrees of freedom. One may equally consider the Dell system 
with arbitrary many degrees of freedom, see below. 
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h{p,q) = f{Q), 
H{P,Q)^F{q). 
Here the functions /(Q) and F{q) are such that 



(7.14) 



dPAdQ^-dpAdq, (7.15) 
which expresses the fact we have an (anti-)canonical transformation. This relation entails that 



dp 



dP 



(7.16) 



At this stage the functions f{Q) and F(q) are arbitrary. However, when the Hamiltonians depend 
on a coupling constant and are such that their "free" part can be separated and depends only on 
the momenta,P| the free Hamiltonians provide a natural choice for these functions: F{q) — ha{q) and 
f{Q)=Ho{Q) where 



H{P,Q)\^.^, = Ho{P). 
With such choice the duality equations become 

ho{Q) = h{p,q), 
Hoiq)^H{P,Q), 

^'^(^'^)i7^(.)^.^(Q)-^^(^'^) 



(7.18) 



(7.19) 



dp 



dP 



Free rational, trigonometric and elliptic Hamiltonians are ho{p) — ft.o(p) — coshp and /io(p) 
cn(p|fc) respectively. 

Note that from the main duality relation. 



it follows that 



which together with (7.19) implies: 



dq 
dP 



dP 



Ho{q)^HiP,Q) 

^ 1 dHjP, Q) 
Q " Hl^iq) dP 

1 dKP: g) 

' h',{Q) dp 



(7.20) 



(7.21) 



(7.22) 



When compared with the Hamiltonian equation for the original system, 

dq ^ dhjp, q) 
dt dp 

^■^ Note that this kind of duality relates the weak coupling regime for h{p, q) to the weak coupling regime for H(P, Q). 
For example, in the rational Calogero case 



(7.23) 



h{p,q) 



2 g2 



(7.17) 



H{P.Q) = — + ^ = — 
^ ' 2 Q2 2 



We recall that, on the field theory side, the coupling constant u is related to the mass of adjoint hypermultiplet and thus 
remains unchanged under duality transformations. 
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we see that P — h'f^{Q)t is proportional to the ordinary time-variable t, while ho^Q) — h{p,q) = E 
expresses Q as a function of the energy E. This is a usual feature of classical integrable systems, exploited 
in Seiberg-Witten theory in the SU{2) case the spectral curve q{t) can be described by 



v\^f^,,y,j^E. (7.24) 

where p is expressed through dq/dt and q from the Hamiltonian equation dq/dt = dH/dp. In other 
words, the spectral curve is essentially the solution of the equation of motion of integrable system, where 
the time t plays the role of the spectral parameter and the energy E that of the modulus. 

Let us consider several simple examples. The simplest one is the rational Calogero. In this case, the 



duality transformation connects two identical Hamiltonians, (7.17). Therefore, this system is self-dual. 
Somewhat less trivial example is the trigonometric Calogero-Sutherland model [^6[ It leads to the 
following dual pair 



Hp, 9) = 



p2 _ Q2 

2 ^ sin^ q~ 2 

(7.25) 



H{P, Q) = Jl - ^ cosP = cosg 



One may easily recognize in the second Hamiltonian the rational Ruijsenaars Hamiltonian jS^, Thus, 
the Hamiltonian rational in coordinates and momenta maps onto itself under the duality, while that with 
rational momentum dependence and trigonometric coordinate dependence maps onto the Hamiltonian 
with inverse, rational coordinate and trigonometric momentum dependencies. Therefore, we can really 
see that the pq-duality exchanges types of the coordinate and momentum dependence. This implies 
that the Hamiltonian trigonometrically dependent both on coordinates and momenta is to be self-dual. 
Indeed, one can check that the trigonometric Ruijsenaars-Schneider system is self-dual: 



h{p, q) = Jl — cosp = cos Q 

(7.26) 



H{P, Q) = Jl - — cosP = cosg 



2z/2 

Now we make the next step and introduce into the game elliptic dependencies. 



8 Two-body system dual to Calogero and Ruijsenaars 

We begin with the elliptic Calogero Hamiltonian^ |3^, 

and seek a dual Hamiltonian elliptic in the momentum. Thus ho{p) = ^ and we seek H{P, Q) = Holq) 
such that ^0(9) = cn{q\k). Eqs.(7.20) become 

T^T + i^?(^' = P-cn'(g|fc) = g^^^^. (8.28) 
Upon substituting 

cn'(q|fc) ^ -sn(q|fc)dn(q|fc) = -^y{l-H^){k'^ + k^m), (8.29) 



(this is because sn^q — 1 — cn^q, dn^q ~ fc'^ + k^cn^q, fc'^ + k^ = 1 and cnq = H) we get for ( ^.28 ) 

'dH_y _ f_ 

dP) ~ 

^^The standard definitions of the Jacobi functions we use below can be found in | |3l| . 



^\ =l-{l-H'){k'' + k^H'). (8.30) 
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Now from the first eqn.( g.2g ) can be expressed tlirough Q and sn^{q\k) = 1 — cu^{q\k) = 1 — H'^ as 

"2 21/2 



so tliat 



^ = 1- 



31) 



Therefore H is an elhptic function of P, namely l3q, 21 



.32) 



H{P, Q) = cn(g|fc) = a(Q) • cn P ^ k'^ + k"^ a"^ {Q) 



with 



ka{Q) 



y/k'^ + k^a^{Q) 



.33) 



In the Hmit v'^ = 0, when the interaction is switched off, a{q) 
cn(P|fc), as assumed in 



.34) 



1 and H{P, Q) reduces to Ho{P) 



We have therefore obtained a dual formulation of the elliptic Calogero model (in the simplest SU (2) 
case). At first glance our dual Hamiltonian looks somewhat unusual. In particular, the relevant elliptic 
curve is "dressed" : it is described by an effective modulus 



ka{Q) 



ka{Q) 



y/k'^ + kM^{Q) ^l-fc2(l_a2(g))' (^-35) 

which differs from the "bare" one fc in a Q-dependent way. In fact k^ff is nothing but the modulus of 
the "reduced" Calogero spectral curve pjt]. 



Now one can rewrite (8.33) in many different forms, in particular, in terms of ^-functions, in hyper- 



elliptic parametrization etc. This latter is rather essential and we refer the reader for the details to the 
paper pl[ |. In particular, this is the hyperelliptic parametrization, where one can conveniently check that 
the symplectic structure dP A dQ is actually equal to dp'''^'^ A da, i.e. to the action-angle variables in the 
Calogero mod el. Here p'^"'^, ( 4.15 ) is the angle along the cycle in the Jacobian and a is the action variable 
given by (1.3). 

All of the above formulae are straightforwardly generalised from the Calogero (rational-elliptic) system 
to the Ruijsenaars (trigonometric-elliptic) system. The only difference ensuing is that the q-dependence 
of the dual (elliptic-trigonometric) Hamiltonian is now trigonometric rather than rational pl[|: 



a'iq) 



trig 



(9) = 1 - 



2i/2 



sinh2 q 



(8.36) 



Rather than giving further details we will proceed directly to a consideration of the double-elliptic model. 



9 Double-elliptic two-body system 



In order to get a double-elliptic system one needs to exchange the rational Q-dependence in (B.28) for 
an elliptic one, and so we substitute a^^^lQ) by the obvious elliptic analogue alii{Q) = 1 — 



Moreover, now the elliptic curves for q and Q need not in general be the same, i.e. k ^ k. 



Instead of (8.28) the duality equations now become 



cn{q\k) = H{P,Q\k,k), 
cniQ\k) = H{p,q\k,k), 



(9.37) 



cn'(g|fc) 



dH{P,Q\k,k) 
dP 



cn' (q\k) 



dH{p, q\k, k) 
dp 
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and the natural ansatz for the Hamiltonian (suggested by ( 8.33| )) is 



H 



H 



{p, q\k, k) = a{q\k, k) ■ en {p P{q\k, k) \ 7(g|/c, fc)^ = acn(p/3|7), 
{P,Q\k,k)=a{Q\k,k)-cn(p l3{Q\k,k) | 7(Q|fc,fc)) = 5cn(P/3|7). 



(9.38) 



Substituting these ansatz into ( |9.37 ) and making use of ( ^.29 ), one can arrive, after some calculations 
[|l),atO 



\q\k,k)=a^{q\k) = 1- 



2v' 



sn^{q\k) 



, I3^{q\k,k) = k'^ + Pa^(q\k), ^^{q\k,k) ^ 



with some constant u and finally the double-elliptic duality becomes pl[ 



ea^{q\k) 

(9.39) 



H{P, Q\k, k) = cn{q\k) = a{Q\k)cn P\lk'^ + k^a'^{Q\k) 



H{p, q\k, k) — cn{Q\k) — a{q\k)cT\ py k'^ + k^a'^{q\k) 



ka{Q\k) 



k'^ + k^a^{Q\k) ^ 



ka{q\k) 



fc'2 + k^a^{q\k)^ 



(9.40) 



(9.41) 



We shall now consider various limiting cases arising from these and show that the double-elliptic 
Hamiltonian ( |9.41 ) contains the entire Ruijsenaars-Calogero and Toda family as its limiting cases, as 
desired. (Of course we have restricted ourselves to the SU{2) members of this family.) 

In order to convert the elliptic dependence of the momentum p into the trigonometric one, the corre- 
sponding "bare" modulus k should vanish: fc ^ 0, k'^ = 1 — k^ ^ 1 (while k can be kept finite). Then, 
since cn(a;|fc = 0) = coshx, 



H^^'^p^q) — >a((?)coshp = i7*'(p,g) 



(9.42) 



with the same 



a^{q\k) = 1 



21.2 



(9.43) 



S'n^{q\k) 

Thus we obtain the SU (2) elliptic Ruijsenaars Hamiltonian. The trigonometric and rational Ruijsenaars 
as well as all of the Calogero and Toda systems are obtained through further limiting procedures in the 
standard way. 



The other limit fc ^ (with k finite) gives a{q\k) — > atrig{q) = 1 



21/-' 
cosh q 



and 



Jldellf^p^ g) 



atrigiq) ■ cn ( p\l fc'2 + k'^a'l^igiq) 



katrigiq) 



= H^"{p, q). 



(9.44) 



This is the elliptic-trigonometric model, dual to the conventional elliptic Ruijsenaars (i.e. the trigonometric- 
elliptic) system. In the further limit of small q this degenerates into the elliptic-rational model with 
atrigiq) ctratiq) = 1 ^ ^2", which is dual to the conventional elliptic Calogero (i.e. the rational-elliptic) 
system, analysed above. 

Our approach has been based on choosing appropriate functions f{q) and F{Q) and implementing 
duality. Other choices o f func tions associated with alternative free Hamiltonians may be possible. Instead 
of the duality relations ( |9.37 ) one could consider those based on /io(p) — sn(p|fc) instead of cn(p|fc). With 
this choice one gets somewhat simpler expressions for (3s and 7^: 

^^For ease of expression hereafter we suppress the dependence of a, /3, 7 on fc and k in what follows using Q-(q) for a{q\k, k) 
and Oi{Q) for a{Q\k, k) etc. 
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A = 1, ls{q\k,k) = kasiq\k), ^^(glfc) = 1 



cn^ {q\k) 



(9.45) 



and the final Hamihonian is now 



Hs{p, q\k, k) = as{q\k) ■ sn(p|fcas((j|fc)). 



(9.46) 



Although this Hamiltonian is somewhat simpler than our earlier choice, the limits involved in obtaining 
the Ruijsenaars-Calogero-Toda reductions are somewhat more involved, and that is why we chose to 



present the Hamiltonian (9.41) first. 

One might further try other elliptic functions for ho{p). Every solution we have obtained by making 
a different ansatz has been related to our solution ( j.41| ) via modular transformations of the four moduli 
fc, k, keff = 7 and kgff = 7. 



10 ^g-duality, r-functions and Hamiltonian reduction 

Now one should extend the results of the two-particle case to the generic number of particles. Unfortu- 
nately, it is still unknown how to get very manifest formulas. Here we describe a method that ultimately 
allows one to construct commuting Hamiltonians of the A^-particle Dell system, but very explicit expres- 
sions will be obtained only for first several terms in a perturbation theory. 

However, in order to explain the idea, we need some preliminary work. In fact, we start from the old 
observation j4l| ^ that, if one study the elliptic (trigonometric, rational) solution to the KP hierarchy, 
the dynamics (dependence on higher times of the hierarchy) of zeroes of the corresponding r-function 
w.r.t. the first time is governed by the elliptic (trigonometric, rational) Calogero system. Moreover, if 
one starts with the two-dimensional Toda lattice hierarchy and studies the zeroes of the r-function w.r.t. 
the zero (discrete) time, the corresponding system governing zeroes is the Ruijsenaars system 

Let us illustrate this with the simplest example of the trigonometric solution of the KP hierarchy. 
This is, in fact, the A^-solitonic solution, the corresponding r-function being 

r({ife}| {i^^,fi,,X^}) = detATxAT (s.j - L,yeS^*'=^'^'-"''^) , L,^ = (10.47) 

Here Vi, jii and Xi are the soliton parameters. Let us impose the constraint Ui = /i^ + e. Then, the 
standard Hamiltonian structure of the hierarchy on the soliton solution gives rise to the Poisson brackets 
(in ^ similar calculation was done for the sine-Gordon case) 



{Xi , Xj } 



{H, - - - e){fii - + e) ■ 



(10.48) 



Note that, upon identification /i; ^ qi, Xi = e^' Y\i 



7^ I ' 



the matrix Lij in (10.47) becomes 



the rational Ruijsenaars Lax operator, with the proper symplectic structure {pi, qj} — Si 
At the same time, (10.47) can be rewritten in the form 



(10.49) 

where x = ti. This determinant is the generating function of the rational Ruijsenaars Hamiltonians. On 
the other hand, the N zeroes of the determinant as the function of ex are just logarithms of the eigenvalues 
of the Lax operator. These zeroes are governed, as we know from ^ by the trigonometric Calogero 
system, and are nothing but the N particle coordinates in the Calogero system: 



N 



T{{tk}\ W^,^l^,X,}) ~ ^e'^^-iJ^ ^ nsin(ea; - q, 



(10.50) 



Therefore, the eigenvalues of the (rational Ruijsenaars) Lax operator are the exponentials e''' of coordi- 
nates in the (trigonometric Calogero) dual system, while the r-function simultaneously is the generating 
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function of the Hamihonians in one integrable systems and the function of coordinates in the dual sys- 
tem. This is exactly the form of duality, as it was first realized by S.Ruijsenaars B. Its relation with 



r-functions was first observed by D.Bernard and O.Babelon (see also [46 



Similarly, one can consider the iV-solitonic solution in the two-dimensional Toda lattice hierarchy. 
Then, 



(10.51) 

This time one should make a reduction Vi = e'lii. Then, Lij becomes the Lax operator of the trigonomet- 
ric Ruijsenaars, and the zeroes w.r.t. to en are governed by the same, trigonometric Ruijsenaars system. 
This is another check that this system is really self-dual. 

Thus constructed duality can be also interpreted in terms of the Hamiltonian reduction |3^, |^ psf . 
In doing so, one starts with the moment map which typically is a constraint for two matrices or two 
matrix-valued functions. In order to solve it, one diagonalizes one of the matrices, its diagonal elements 
ultimately being functions (exponentials) of the coordinates, while the other matrix gives the Lax operator 
(its traces are the Hamiltonians) . Now what one needs in order to construct the dual system is just to 
diagonalize the second matrix in order to get coordinates in the dual system, while the first matrix will 
provide one with the new Lax operator (its traces are the Hamiltonians in the dual system). 

As an illustrative example, let us consider the trigonometric Calogero-Sutherland system |^7|] . In 
order to get it, one can start from the free Hamiltonian system H = tiA^ given on the phase space {A, B) 
of two N X N matrices, with the symplectic form tr(i5A A B~^5B (i.e. A lies in the gl{N) algebra, while 
B lies in the GL{N) group). Then, one can impose the constraint 

A - BAB-^ = v{l - P) (10.52) 

where i/ is a constant, I is the unit matrix and P is a matrix of rank 1. Now one has to make the 
Hamiltonian reduction with this constraint. First of all, one solves the constraint diagonalizing matrix 
B and using the gauge freedom of the system to transform P into the matrix with all entries unit. It 
gives|^ Aij = Pi6ij + (1 — Sij) , where bi are the diagonal elements of B. This is nothing but the Lax 
operator of the trigonometric Calogero-Sutherland model. Traces of powers of this Lax operators give 
the Calogero Hamiltonians. 



One can also diagonalize the other matrix, A. Then, (10.52) can be rewritten in the form 



AB - BA = u{B - P) (10.53) 

where P is another rank 1 matrix which is can be gauged out to the matrix with all entries depending only 
on the number of row (let denote them Xi/v). Solving then this constraint, one obtains B — ^ 
where are the diagonal elements of the matrix A. This is nothing but the Lax operator of the rational 
Ruijsenaars-Schneider system, traces of its powers give a set of the dual Hamiltonians. 

Certainly, this picture nicely explaining the Ruijsenaars observation of duality is just the r-function 
approach told in different words. 



11 Dual Hamiltonians for many-body systems 

Unfortunately, the scheme described in the previous section does not work for elliptic models with arbi- 
trary number of particles. The reason is that even having the dual Lax operator calculated, one is typically 
not able in elliptic cases to construct the dual Hamiltonians, i.e. invariant combinations involving the 
Lax operator. 

Instead, in order to obtain dual Hamiltonians, here we use the other approach discussed above, that 
dealing with zeroes of the r-functions [|l], This approach gives no very explicit form of the dual 
Hamiltonians. Therefore, we construct later a kind of perturbative procedure for the Hamiltonians which 
allows one to get them absolutely explicitly term by term [E2[ . 



Since (10.52) does not fix the diagonal elements of A, they are parametrized by arbitrary numbers pi 
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Consider SU{N) {N — g + 1) system. The whole construction is based on the fact that the spectral 
curve of the original integrable system (Toda chain, Calogero, Ruijsenaars or the most interesting double 
elliptic system) has a period matrix Tij{a), i, j — I, . . . , N with the special property: 



N 

^T,j(a) = r, Vi (11.54) 



where r does not depend on a. As a corollary, the genus- A'^ theta-function is naturally decomposed into 
a linear combination of genus-.g thcta-functions: 



TV N \ N-1 



e(^)(K|r„) = ^ exp ^ r„n,n, + 2^:^Y,n^P^ = J2 ^[-^.o]i^C\NT) ■ Qi'\p,\f,j) (11.55) 

MiGZ y i,j=l j=l / A:=0 

where 

(N ^ \ 

-ITT 4 + (11.56) 

= l 1=1 J 

and p, = (+p^, Pi = 0; = + j^, J2f=i ^ij = 0, Vi. 

Now we again use the argument that zeroes of the KP ( Toda) r-function (i.e. essentially the Rie- 
mannian theta-function), associated with the spectral curve ( |13.73D are nothing but the coordinates q. 



of the original (Calogero, Ruijsenaars) integrable system. In more detail, due to the property ( 11.54 ), 
0(^)(p|7^) as a function of C = 12iLi Pi elliptic function on the torus (1, t) and, therefore, can be 
decomposed into an A^-fold product of the genus-one theta-functions. Remarkably, their arguments are 
just C - 9*: 



N 



e(^)(p|T) = cip, T, t) n 0(C - (7.(P, T)\t) (11.57) 
1=1 

(In the case of the Toda chain when r — > ioo this "sum rule" is implied by the standard expression for 
the individual e'* through the KP r-function.) Since one can prove that qi form a Poisson-commuting 



set of variables with respect to the symplectic structure (11.58) |Q, this observation indirectly justifies 
the claim |2|] that all the ratios Qk/^i (in order to cancel the non-elliptic factor c(p, T, r)) are 
Poisson-commuting with respect to the Seiberg-Witten symplectic structure 



N 



^dpi^^hdai (11.58) 



i=l 



(where = —ai — ... — oat^i, i.e. J2i=i ^ 

^,^1=0 VA;,;,m,n (11.59) 

The Hamiltonians of the dual integrable system can be chosen in the form Hk = 8fc/8o, k = 1, . . . ,g. 
However, these dual Hamiltonians are not quite manifest, since they depend on the period matrix of 
the original system expressed in terms of its action variables. Nevertheless, what is important, one can 
work with the 0- functions as with series and, using the instantonic expansion for the period matrix Tij, 
construct Hamiltonians term by term in the series. 



^*It can be equivalently written as 

0^{^J,ek} + e,{^k,e^} + ek{e,.,eJ} = \/i,j,k 
{loge^ioge^} = <j iog|^,iog0fc \ vi,j,k 
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12 Dual Hamiltonians perturbatively 



Let us see how it really works. We start with the simplest case of the Seiberg-Witten family, the periodic 
Toda chain. 



12.1 Perturbative approximation 



First we consider the pertm'bative approximation which is nothing but the open Toda chain, see (3 
and (6.9). This is the case of the perturbative Ad pure N — 2 SYM theory with the prepotential 



1 ^ 



(12.60) 



i<3 



Then, the period matrix is singular and only finite number of terms survives in the series for the theta- 
f unction: 



e(^)(p|r) 



JV-l 
fe=0 



.27rifcC TjiO) 



Hf\p,a) 



E 

I,[I\ = k i&I 



Here 



A 



QCD 



(12.61) 



(12.62) 



(12.63) 



and / are all possible partitions of N indices into the sets of fc = [/] and N — k = [I] elements. Pa- 
rameter Aqcd becomes significant only when the system is deformed: either non-perturbatively or to 
more complex system s of the Calogero-Ruijsenaars-double-elliptic family. The corresponding r-function 
q(^^p\T) , eq.( |l2.6l" ), describes an iV-soliton solution to the KP hierarchy and is equal to the determinant 



( [10.47[ ). The Hamiltonians fff \ eq.( |l2.62| ) dual to the open Toda chain, are those of the degenerated 



rational Ruijsenaars system, and they are well-known to Poisson-commute with respect to the relevant 
Seiberg-Witten symplectic structure ( 11.58|) . This is not surprising since the open Toda chain is obtained 
by a degeneration from the trigonometric Calogero system which is dual to the rational Ruijsenaars. 

The same construction for other perturbative Seiberg-Witten systems ends up with the Hamiltonians 
of more sophisticated systems. 

For the trigonometric Calogero system (perturbative 4d A/" = 4 SYM with SUSY softly broken down 
to A/" = 2 by the adjoint mass M) the Poisson-commuting (with respect to the same (11.5S)), whose 



perturbative prepotential is given by the sum of (5.2) and (5.3), Hamiltonians iJ^ 
with 



(0) 



are given by (12.62) 



Af2 



(12.64) 



i.e. are the Hamiltonians of the rational Ruijsenaars system which are really dual to the trigonometric 
Calogero model. 

For the period matrix (effective charges) ( |5.5| ) of the trigonometric Ruijsenaars system (perturbative 
bd M = 2 SY M com pactified on a circle with an e twist as the boundary conditions) the Hamiltonians 
are given by ( 12.62 ) with 



yjsm\i[aij + e) sinh(ajj — e) 
sinh an 



(12.65) 



i.e. are the Hamiltonians of the trigonometric Ruijsenaars system, which is, indeed, self-dual. 

Finally, for the perturbative limit of the most interesting self-dual double-elliptic system (the explicit 
form of its spectral curves, whose period matrix is supposed to be (5.7), is yet unknown but in the case of 
two particles, see the next section) the relevant Hamiltonians are those of the elliptic Ruijsenaars system, 
given by the same ( 12.62| ) with (cf. with (5.7)) 
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where r is the modulus of the second torus associated with the double elliptic system and are a^- 
rescaled, pl[ |. 

In all these examples H^^^ = 1, the theta- functions G*-^' are singular and given by determinant 
(solitonic) formulas with finite number of items (only terms with rii = 0,1 survive in the series expansion 
of the theta- function) , and Poisson-commutativity of arising Hamiltonians is analytically checked within 
the theory of Ruijsenaars integrable systems. 

12.2 Beyond the perturbative limit 

Beyond the perturbative limit, the analytical evaluation of O*-^' becomes less straightforward. 



The non-perturbative deformation (4.13) of the curve (3.8), is associated with somewhat sophisticated 
prepotential of the periodic Toda chain, 

-.AT oo 



i<j k=l 



The period matrix is 



-e 



i=ij k=ii,j i k^i 



+ .. 



with 



■' \ daidaj 



Z,, = = Z^> { 1 - »^ +...), i^j (12.70) 



(12.68) 



and T in this case can be removed by the rescaling of A. Then, 



AT 

e^^**^^'~^^^ Zfj^ZfiZjj^Zji Y\ ^imZjjnZkrnZlm + ■ ■ ■ 1 + 

i^j^k^l m^i.jjkjl j 

/ w (12.69) 
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The first few corrections -E^'^' to the prepotential are exphcitly known in the Toda-chain case [Q, for 
example, 



1 ^ 2 

^^^^=-^En(4°') (12.71) 



2i7r 

i—l k^i 



The coefficients Ok in ( |12.69 ) are expanded into powers of (A/a)^^ and the leading (zeroth-order) 



terms are exactly the perturbative expressions (12.62). Thus, the degenerated Ruijsenaars Hamiltonians 



(12.62) are just the perturbative approximations to the Hk — Qk/^o ~ the full Hamiltonians of the 
integrable system dual to the Toda chain. 

The same expansion can be constructed for the other, elliptic systems. Note that the period matrix 
in these cases is expanded into powers oi q = e^'^*'^: 



fe=l 



a,- 



and the dimensional constant /i depends on the system. E.g., in the Calogero model ^ — iM etc. The 
limit to the Toda system corresponds to g ^ 0, (7/^^^=fixed. 

Note that the first (instanton-gas) correction to the prepotential is indeed known to be of the same 



universal form (12.71) not only for the Toda chain, but also for the Calogero system |j49|. For the 



Ruijsenaars and double-elliptic systems eq.(12.71) is not yet available in the literature 



In p2] the perturbative expansions for the Hamiltonians were tested for Poisson commutativity, up 



to few first terms and small number of particles N . In particular, it was checked that (12.71) is true for 
both the Ruijsenaars and Dell systems. 



The main lesson we could get from this consideration is that if the universal expressions like (12.71) 



in terms of perturbative z[^^ , the same for all the systems, will be found for higher corrections to the 
prepotential^^, this will immediately give an explicit (although not the most appealing) construction of 
the Hamiltonians dual to the Calogero and Ruijsenaars models and, especially important, the self-dual 
double-elliptic system which was explicitly constructed above only for N — 2. 



Let us note now that the SU (2) Hamiltonian (9.41 ) has a rather nice form in terms of the p, q variables 



or their duals P, Q. However, this our general SU{N) construction of the Hamiltonians as ratios of genus 
g 0- functions uses another kind of canonical variables - angle- action variables pf"*^, Oj. The flat moduli 
play the central role in Seiberg-Witten theory theory, while Qi are rather generalizations of the algebraic 
moduli. 

It is an interesting open problem to express the Hamiltonians in terms of Pi and Qi, perhaps, they 
can acquire a more transparent form, like it happens for SU{2). Another option is to switch to the 
"separated" variables Vi, Qi such that §^_VidQi — aiSij (while generically — §j^.^iPidQi). 
Generically, they are different from Pi, Qi but coincide with P, Q in the case of SU{2). 



13 Dell systems: 6d adjoint matter 



Thus, we have constructed the Dell system and a set of dualities that acts on the Toda-Calogero- 
Ruijsenaars-Dell family. They are all put together in Fig. 2. 

To avoid possible confusion, the recurrent relation | po| , [l4| for the Toda-chain prepotential, 

a^jf 92^ 92 



91ogA2 dlogAdai 91ogA9aj dpidp 



■logeo 



p=0 



does not immediately provide such universal expressions. Already for Z^^"i this relation gives 



N 

E 

r(0) 



which coincides with (12.71) for the Toda-chain -Z^j , eq.{12.63), but is not true (in variance with (12.71)) for Calogero 



AO) 



eq.(12.64). Meanwhile, the recurrent relations of rcf. for the Calogero system provide more promising expansion. 
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Figure 2: Action of the coordinate-momentum duality on the Calogero-Ruijsenaars-Dell family. Hooked 
arrows mark self-dual systems. The duality leaves the coupling constant intact. 



Since the theories of this family describes low energy limits of different SUSY gauge theories, this 
table describes the dualities between different gauge theories too. Say, perturbative limit of the M 
theory (trigonometric Calogero) is dual to a special degeneration Q of the perturbative limit of the 5d 
theory (rational Ruijsenaars), while full, non-perturbative 5d theory (elliptic Ruijsenaars) is dual to the 
perturbative limit of the 6d theory (Dell system). 

To conclude our discussion of dualities and Dell systems, we calculate the perturbative prepotential 
of the 2-particle Dell system ( ^.41 ) and show that it coincides with what is to be expected for the 6d 



theory with adjoint matter, supporting the identification of the Dell system with the 6d gauge theory. 
Let us start from the full Dell model for 2-particles. Then, it is given by the Hamiltonian ( p. 41 ) 



parametrized by two independent (momentum and coordinate) elliptic curves with elliptic moduli k and 
k. As usual, for two particles we can construct the full spectral curve from the Hamiltonian (H, see also 



section 4 and (7.24)), 



iJ(C, ^\k,k) = u ( = cn(g|fc) ) (13.73) 
It is characterized by the effective elliptic moduli 

ka{q\k) J ka{q\k) n-i7A\ 

keff — —, keff — (13.74) 

I3{q\k,k) f}{q\k,k) 



where the functions a and (3 are manifestly given in ( 9.39| ). Coordinate-momentum duality interchanges 
k k, keff keff (and q,p ^ Q,P), while in general SU{N) case, the model describes an interplay 
between the four tori: the two bare elliptic curves and two effective Jacobians of complex dimension 
g = N-l. 

The generating differential in 6d theories is of the form dS — (d£,, with ^ living on the coordinate 
torus and C ~ on the momentum torus. 

Now we calculate the leading order of the prepotential expansion in powers of k when the bare spectral 
torus degenerates into sphere. In the forthcoming calculation we closely follow the line of pO|] . 

When fc 0, S'n{q\k) degenerates into the ordinary sine. For furth er convenience, we shall parametrize 
the coupling constant 21/^ = sn^(e|fc). Now the spectral curve (13.73) acquires the form 



sn2(e|fc) 



sin^ ^ 



cn{C(3 



-eff 



(13.75) 



Here the variable ^ lives in the cylinder produced after degenerating the bare coordinate torus. So does 
the variable x = 1/ sin^ ^. Note that the A-period of the dressed torus shrinks on the sphere to a contour 
around a; = 0. Similarly, i?-period can be taken as a contour passing from a; = to a; = 1 and back. 

The next step is to calculate variation of the generating differential dS = rjdS, w.r.t. the modulus u in 
order to obtain a holomorphic differential: 
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where f/^ = ^.Vx ■ Since 



9a _ f ddS _ r 1 

we deduce that m = cn(a|A:) and [/ = . The ratio of the B- and A-periods of dv gives the period 

matrix 



U f dx 1 , / /t\ 1 , C/2 , ^ 

T = — = = -— hm log- + —log- — (13.78 

Trio xV(a;-l)(f/2-x) «7r«^oV ^4/ nr "^l-U^ ^ ' 

where k is a small- a; cut-off. The U dependent part of this integral is finite and can be considered as the 
"true" perturbative correction, while the divergent part just renormalizes the bare "classical" coupling 
constant r, i.e. classical part of the prepotential (see for further details). Therefore, the perturbative 
period matrix is finally 



' finite 



i sn^(e|fc) — sn^(a|fc) 



-log 

TT 



, const — > -log ^271 (13.79) 

sn^(a|A;) tt "iW 

and the perturbative prepotential is the elliptic tri-logarithm. Remarkably, it lives on the bare momentum 
torus, while the modulus of the perturbative curve (13.75) is the dressed one. 



v{a-e) 
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